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Abstract. 

The glow discharge's main ionization breakdown processes have been understood 
for about one hundred years. The glow discharge, however, still remains as an avenue of 
fresh research in pattern formation and far from equilibrium systems. The primary and 
secondary ionization processes can be mathematically modeled as general branching 
processes. Not only is the Townsend breakdown criterion obtained but the ionization 
avalanche can be characterized as a branching set with a unique Hausdorff fractal 
dimension. These fractal dimensions show applicability using similarity principles and 
Paschen's Law. 
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The cylindrical geometry, planar electrode DC glow discharge is one of the most 
familiar features in the study of ionized gases. Though the features of the DC 
discharge have been described since at least the time of Michael Faraday, they were 
not theoretically understood until the work of John S.E. Townsend early in the 20th 
century [TJ E]. Townsend was able to describe the gas breakdown in terms to two 
ionization cascades. A primary ionization cascade of electrons in the gas as voltage is 
increased due to electron avalanches beginning with electrons from atoms ionized by 
external radiation like UV sources, natural radioactivity, or cosmic rays. A secondary 
ionization cascade is caused by electron ionizations from the positive ion counterparts 
to the electrons colliding with other atoms in the cathode or other gas atoms. At 
breakdown, the discharge becomes self-sustaining, transitioning from the dark Townsend 
discharge to the luminous glow discharge. At this point secondary ionization provides 
enough initial electrons and external radiation is no longer needed. This has been long 
understood, however, the glow discharge has recently become a fertile area of research 
into pattern formation and non-equilibrium systems [31 SE]- In particular, Sijacic and 
Ebert [HI E] have formulated the Townsend breakdown in terms of a supercritical and 
subcritical bifurcation. This paper will investigate the ionization cascades in terms of the 
topology of the ionization cascade. Specifically, by formulating the ionization avalanche 
as a branching process, it can be described with a Hausdorff fractal dimension. 

1. General Branching Process 

The theory of branching processes is a theory of the growth of a population through 
generations given each member of the population leaves the population, remains single, 
or produces a number of progeny given the transition probabilities of a fixed Markov 
process. It was first described by Francis Galton, H.W. Watson, and earlier by 
Bienayme. Though related research in the probabilistic nature of breakdown has been 
done [HI E] , in this paper we will use the general branching process described by Harris 
[TT] to derive the Townsend breakdown condition. 

The objects considered by our branching process will be the electrons in the 
gas discharge primary and secondary ionizations. The transition probability p& is 
the probability that an electron will produce k electrons through collision from one 
generation to the next. The probability generating function f(s) is 



Where /'(l) is the defined as the expectation value, m, which is the average number 
of new electrons each individual electron will produce each generation. If m > 1 there 
is a non-zero probability of an infinite number of electrons in the limit of the number 
of generations goes to infinity. According to the theory of branching processes, as the 
number of generations goes to infinity, the electron population trends to either infinity 
or zero. 



oo 
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The general branching process is a multitype branching process where an attribute 
of the members of the population are described by using a random finite set of variables, 
x, which belongs to an interval X, a subset of the positive real numbers. This variable 
can be used to describe among other attributes, the age, energy, or in this paper, the 
position of a particle. The point distribution uj is the set of all distinct points x\ to x k 
with positive integer weights n\ to n k . 

uj = (x 1 ,n 1 ;x 2 ,n 2 ; ...x k ,n k ) (2) 

The integer weight n k represents the number of objects point x k . Define Qx as the set 
of all point distributions. Define a set function uj where £o(A) = J2 x eA n i- If A is null 
oj(A) = 0. In addition each subset A can have a probability defined P(u(A)) = r which 
is the probability that there are r total objects in the subset A. Finally, a function s(x) 
is defined 

s(x) — / s(x)du(x) — nis(xi) +n 2 s(x 2 ) + ...n k s(x k ) (3) 

J X 

Where s(x) is 1 if x is in A and otherwise. The moment generating functional 
(MGF), the function that supplements the probability generating function in point set 
distributions, is defined as 

<S>( s )=Ee-f sdw =[ e~f sdw dP (4) 
Jn x 

but also has an alternative formulation 

=/o + E /•••/ fn(*i, ...,x n )e- s ^-- s ^dx 1 ...dx n (5) 

n=l 

where f is the probability of no objects and f n (xi, ...,x n )dx\...dx n is the probability of 
n objects in which the first is in the interval {x\,x\ + dx\) etc. This is the formulation 
we will use in our derivation. The purpose of defining the MGF is to enable us to use 
each x k and its corresponding weight as the number of electrons at a point x k on the 
length of the ionization region from to d, the length of the cathode sheath. For a 
general branching process the corresponding variable to the expectation value, to, is the 
eigenvalue p of the first moment of the MGF where the first moment is defined as 

M 1 (x,A) = E x Z 1 (A)= ~ 9 *£ 8A) (6) 

at t=o 

Zi(A) is the number of objects in A in the first generation. The first moment also has 
a density m\ defined as 

M 1 (x,A) = [ mi (x,y)dy (7) 

J A 

The density mi is also related to the eigenvalue, p which plays the same role in 
multitype branching processes as the expectation value, m, and the right eigenfunction 
fj,(x) and the left eigenfunction u(y) by 

m 1 (x,y) = p^x)u(y)[l + 0(A)] (8) 
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1.1. Primary Ionization 

For the primary ionization in the Townsend discharge we define the transition 
probabilities of po being the probability of recombination, pi the probability of 
scattering, and p<i the probability of ionization. An electron is assumed to only be 
able to ionize a maximum of one other electron in a collision. 

Let A be the mean free path for electrons in the gas. Therefore e~* be the 
probability of no collision of an electron over the distance d jTHj. Here we consider 
d to be the length of the ionization region in the cathode sheath. The probability of an 
electron starting at point x colliding at a distance y is e ~ . We can then define the 
MGF for the first generation as 

(d-x) f (y-x) 



= (e * ) + £p» / e-^e-^dy (9) 

n=0 JX 

which is expanded as 

ti\ (d-x) I'd (y-x) „ fd (y-x) , I'd (y-x) „. 

$W(s) = (e ~)+p / e ~e°dy+p 1 / e ~e' tSA dy+p 2 / e ~e~ 2tSA dy(10) 

Jo Jo Jo 

sa, given the equation above, is the total number of electrons in the first generation. 

For the primary cascade we assume one initiating electron so — 1. The first moment 

as defined in El is 

fd (y-x) fd (y-x) 

M 1 =p 1 e —dy + 2p 2 e ~dy (11) 
Jo Jo 

which can be simplified as 



d 



(y-x) 



M 1 = m e ~dy (12) 
Jo 

where m if the expectation value of the probability generating function of the primary 
ionization cascade. 

The density is then 

(y — x) 

mi = me ~ (13) 

Finally, given |H] we can see the values of the eigenfunctions must be n{x) = e^ and 
v{y) = e~A so the eigenvalue, and thus expectation value, of the general branching 
process is 

p = m (14) 

Thus the eigenvalue p equals f'(l) = m of the transition probabilities as expected. 
The expected value of the number of objects, Z, in generation n is given by Z n = m n 

According to Townsend's ionization theory, the number of electrons at a length d 
created from one primary electron is = e ad — 1 or iV d ~ e ad . However, since our 
branching process assumes a maximum of two electrons produced every generation, the 
correct base for the exponential growth is 2 not e. If d is again the total length of the 
ionization region in the cathode sheath, 

m n = 2 ad (15) 
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1.2. Secondary Emission 

Though secondary emission is initiated by ions, since there is an equality of the number 
of ions and primary electrons we can assume the primary electrons to be the causal 
objects in the secondary emission. The transition probabilities for secondary emission 
are slightly different. Pq is the probability an ion (primary electron) collides with the 
cathode and is absorbed. P\ is the probability an ion (primary electron) collides with 
the cathode, is absorbed and emits a secondary electron. P\ is equal to 7, Townsend's 
secondary ionization coefficient. 
Then MGF is then 

$« = Po [ d e °dy + Pl [ d e- ts{y) dy (16) 
Jo Jo 

and the first moment and density are 

M 1=Pl s A [ d e°dy (17) 
Jo 

mi = msA (18) 



However, in this case sa is not one since we are not assuming a single parent, rather 
we are assuming all the ions from the primary electrons at all points to d to be initial 
objects in secondary emission so sa = Z n . Given |H] we can see that neither x nor y are 
in the first moment's density so the left and right eigenf unctions must both be equal to 
1 giving us a eigenvalue p 2 of 

p 2 = ms A (19) 

Given we showed earlier Z n = 2 ad — 1 and m = 7 we see 

p 2 = 7(2 ad - 1) (20) 

Since an infinite avalanche only happens when nonzero probability if p 2 > 1 we see 
that we have the Townsend breakdown condition under conditions of binary division. 



2. Branching Trees 

Now that we have shown the Townsend breakdown criterion to be a result of a general 
branching process, we can extend this result to branching trees. Let a branching tree, 
T, be defined as a graph G where the set of vertices V(G) consists of all objects in all 
generations and let the set of edges E(G) be the set of directed edges point from an 
object to its immediate "descendants" in the next generation. The branching tree starts 
with a single root at generation and expands with the reproduction of the objects in 
subsequent generations. 

There has been much research linking probabilistic branching processes and trees 
[T2"l [T5J Ej They have conclusively shown that a branching process can be represented 
by a branching tree. 
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Figure 1. A branching tree with four generations starting from its root. 

The branching tree has a boundary, dT, which is the set of all infinite rays from 
the root. In other words, the set of all branches on the branching tree which do not 
include leaves. dT forms a space on which lie the objects in the branching process. The 
distance metric is defined for two points x, y in dT as e~ xAy where x A y is the index of 
the closest generation in which x and y have a common ancestor. An interesting result 
is given this definition of a metric space, the branching tree forms a random fractal with 
a fractal dimension defined as 

D = log br[T] (21) 

6r[T], the branching number, is defined in general as the supremum of a variable 
A where A - " defines the capacity for flow at an edge n generations from the root. If A 
holds a value greater than br[T] there can be no flow in the tree. For trees generated 
by branching processes, br[T] is equal to the expectation value of the branching process 
(p). One consideration before extending the branching tree results to the primary and 
secondary ionization of the glow discharge is the fact that the theoretical branching 
tree model is a set of infinite rays while the breakdown processes in the glow discharge 
occur only over a finite length. It is my contention that since the ionization cascade is 
limited by an external physical constraint, the current limiting resistor in series with 
the discharge tube, that an infinite theoretical model can provide a close approximation 
for the ionization cascades which are supercritical as both the primary and secondary 
ionization are at breakdown. 

Given the above equation for fractal dimension we can find expressions for the 
fractal dimension of both the primary and secondary ionization avalanches. For primary 
ionization we first use the expression We then solve for log m to get 

D 1 = aAln2 (22) 

Where A is the mean free path of the gas. We assume the mean free path is the average 
distance an electron travels before a collision and is thus 



A = d/n 



(23) 
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Where n is the number of generations in the branching process. Also you can interpret 
Di as 

Di = (A/Aj) In 2 (24) 

Where A, is the average ionization length. Di will subsequently be known as the primary 
fractal dimension. For secondary ionization since we know m is the Townsend breakdown 
condition, we obtain 

D 2 = log 7 + ad In 2 (25) 

D 2 can also be expressed in terms of Di 

D 2 =\og 1 + nD l (26) 

D 2 will subsequently be known as the secondary fractal dimension. Obviously D 2 
scales with the size of the system while Di is dependent only on the characteristics of 
the gas and electric fields that give rise to the first ionization coefficient. 



3. Similarity Principles and Paschen's Law 

The similarity principles are a common concept in dealing with ionized gases. Similarity 
principles allow us to compare discharges of different linear dimensions. With similarity 
relations we can relate various macroscopic variables of the two discharges with the 
same potential difference drop to each other. The fractal dimensions of the primary and 
secondary ionization cascades both yield surprising results given the similarity principles 
Assume two discharges. One has a length d\ while the other has a length d 2 where 
d 2 = kd\. Assume both have the same cylindrical geometry and planar electrodes and 
that the voltages between the electrodes in each are U\ = U 2 . Two well known similarity 
principles ^3] are the similarity principle for the mean free path, A1/A2 = 1/k and the 
first ionization coefficient ai = ka 2 . We can see that there is a corresponding similarity 
relation for the fractal dimension of the first ionization cascade 

aiAi = a 2 \ 2 (27) 

Therefore, if two discharges are similar, they have the same primary fractal 
dimension and thus share an identical underlying topology. Therefore, with two 
discharges with the same fractal dimension, all other similarity principles necessarily 
follow including p\d\ = p 2 d 2 . So two discharges with identical pd and U also have the 
same primary fractal dimension. Thus there exists a relationship Di = f(pd). Given 

a = Ape- Bp / E (28) 



And given the mean free path equation 
. kT 
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Where k is Boltzmann's constant, T is the temperature of the gas and a is the 
cross-section. Assuming E = U/d we can deduce f(pd) of being of the form 

MT -Bpd/U 

D x = (30) 

a 

This result explains in greater clarity the relationship between the primary fractal 
dimension and the similarity principles. Given any discharge with the same gas 
constants, the value of pd/U is the same for all discharges with the same primary 
fractal dimension. Thus the similarity principles of p\d\ = p 2 d 2 given U\ = U 2 are a 
subset of a larger similarity principle. The principle that all discharges with the same 
primary fractal dimension follow 
Pidi = p 2 d 2 
Ui ~ U 2 

3.1. Paschen's Law 



(31) 



Paschen's Law is well-known as the principle that the breakdown voltage of a gas is 
constant for a constant pd. All points on the Paschen curve indicate the voltage where 
the gas the curve describes (assuming planar electrode geometries) transitions from the 
Townsend discharge to the normal glow. Given the results in|^we can now determine 
that for any discharge with a given pd, its breakdown voltage occurs at the same primary 
fractal dimension. In addition, since = C\ where C is a constant and thus U = C 2 pd 
where C\ = 1/C 2 we can deduce that on any Paschen curve plot by drawing a straight 
line of any positive slope from the origin, all points the straight line intersects on the 
Paschen curve have the same primary fractal dimension on breakdown. 

The secondary fractal dimension also provides a new insight into Paschen's Law. 
By definition at every point on the Paschen curve the Townsend breakdown criterion is 
satisfied. In fact, theoretical expressions for the Paschen Curve are derived using the 
breakdown criterion as a starting point. Loeb ^0] in fact describes the condition for a 
stable discharge according to Holst-Oosterhuis-Seeliger as 

(1 + -) = e^ adx (32) 
7 

Once the normal glow is initiated, if the gas has a value of pd greater than or 
equal to the value of pd at the Paschen minimum, the gas voltage and pd values adjust 
themselves to the values at the Paschen minimum. If the value of pd was less than 
the minimum, the gas maintains itself at the higher breakdown value of voltage and the 
lower pd value. These conditions only change once voltage and current return to positive 
differential conductivity in the abnormal glow. Thus, given that the normal glow always 
maintains itself on the Paschen curve under typical conditions the normal glow can be 
described in all states by the secondary ionization fractal dimension D 2 = 0. This 
reinterprets Paschen's Law as a contour line representing D 2 = 0. This is significant in 
that it shows the normal and abnormal (due to the stability condition) glow despite its 
voltage or pd product, is uniquely defined by the underlying topology of its secondary 
ionization. 
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4. Discussion 

A key question raised by the foregoing is whether fractal dimension in the discharge 
ionization cascades has any additional meaning besides a restatement of the similarity 
principles and the Townsend breakdown criterion. The author would argue that it opens 
a new avenue for investigation by helping us understand that the similarities between 
discharges due to their primary and secondary fractal dimensions helps us to more 
clearly understand or interpret the underlying ionization processes in discharges that 
would seemingly have no similarity, even using the standard similarity principles. The 
fractal dimensions discussed can be easily calculated given experimental data for the 
needed variables. Future suggested research would include experimental measurements 
of both fractal dimensions under different conditions to ascertain whether phenomena 
are similar to or dissimilar from those postulated in this paper. One key point the 
author would like to raise is the scaling that leads to the fractal dimensions mentioned 
in this paper is not related to the scaling in avalanche size in critical phenomena. 
Since it deals with theoretically infinite avalanches it does not answer nor ask the 
same questions of avalanche phenomena as theories such as self-organized criticality. 
It would be interesting, however, to see these techniques applied to a wider class of 
critical phenomena. I would like to thank Robbie Stewart and Steve Hansen for help 
and suggestions regarding experimental apparatus and Russell Lyons and Yuval Peres 
for helpful comments on the mathematics of branching sets. 
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